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1 Introduction and statement of the main result 

Let M be a smooth orientable manifold of dimension n endowed with a volume form 0. We 
denote by C°°(M) and X°°(M) the spaces of smooth real- valued functions and vector fields 
on M respectively. We shall only work in the smooth (C°^) category, although most of our 
results can be easily adapted to lower regularity. 

This letter is focused on the integrability properties of vector fields. More precisely, we 
are interested in proving some sufficient conditions which imply the existence of explicit first 
integrals (i.e. computable). We recall that a function H S C°^{U), where [/ C M is an open 
subset of M, is a first integral of the vector field X in ?7 if Lx{H) = 0, which implies that H is 
constant along the integral curves of X. As usual, Lx denotes the Lie derivative with respect 
to the vector field X. We recall that an equivalent way of writing the first integral condition is 
Lx{H) = ix{dH) = X[H) = 0, where d is the exterior derivative and i denotes de contraction 
operator. In this paper we are interested in global or semi-global first integrals in the sense that 
the open set [/ C M where the first integral is well defined is known a priori. 

The importance of the existence of a non-constant first integral lies in the fact that the 
trajectories of the vector field X leave invariant the level sets of the function H, and hence this 
is a strong constraint on the dynamical behavior of the vector field (e.g. it prevents from the 
existence of ergodic trajectories on M). Unfortunately, it is generally very difficult to compute 
a first integral (whose existence is, in fact, a rather non-generic phenomenon). This is the reason 
why it is necessary to introduce some auxiliary objects to study the integrability properties of 
a vector field. 
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There are two important tools which can be used to analyze the existence of first integrals: 
integrating factors and normalizers. We say that a function / E C°°(M) is an integrating 
factor of the vector field X (with respect to a volume form 0) if divQ{fX) = 0. Recall that 
the divergence of a vector field X with respect to a volume form J7, divQ(X), is defined by 
Lx{^) = diYn{X)Q. A vector field Y G je°°(M) is a normalizer of X if [X,Y] = XX for some 
function A G C°°{M) (the notion of normalizer was introduced by Lie and Engel, see e.g. the 
recent translation [14, Theorem 20, Chapter 8]). Remember that [•, •] stands for the Lie bracket, 
which is defined as [X, Y] = LxiY). Normalizers are also called orbital symmetries, and in the 
case that [X, Y] = 0, it is usually said that the vector field y is a centralizer (or symmetry) 
of X. 

There is an abundant literature on the connection between first integrals, integrating factors 
and normalizers. The fact that the existence of normalizers gives rise to computable first 
integrals in terms of quadratures is well known since Lie's works, see e.g. the monographs 
on the subject by Olver [16] and Bluman-Anco [3] or the comprehensive geometric account of 
reduction and integrability by Sherring and Prince [20]. Some recent results on the relationship 
between symmetries and first integrals can be consulted in [4, 13, 17, 19] and references therein. 
On the other hand, the existence of integrating factors is not generally enough to ensure the 
existence of first integrals, except for the case that M has dimension 2 [8, Chapter 2.11]. In the 
case that M = and the vector field X is polynomial, some algebraic algorithms have been 
proposed to compute polynomial integrating factors and elementary first integrals [18, 21]. For 
volume-preserving vector fields, i.e. divQ(X) = (a particular case of vector fields admitting 
an integrating factor) which have volume-preserving normalizers, there are several procedures 
to compute explicit first integrals using quadratures [9, 11, 15]. We also note that integrating 
factors are also directly related to normalizers, and in fact can be obtained from them under 
suitable hypotheses, see e.g. [5]. 

In this letter we follow a totally different approach, which combines the existence of inte- 
grating factors with a suitable extension of the notion of normalizer (see (z) below). The main 
difference with the approaches explained in the previous paragraph is that we provide an explicit 
formula for the first integral of the vector field in terms of the normalizer and the integrating 
factor, and there is no need of quadratures, reduction of order or algebraic manipulations to 
get this first integral. Previous results in this direction were obtained in [6, 7, 10] for volume 
preserving vector fields with centralizers (see Section 3 for more details). Let us now state the 
main theorem of this letter. 

Theorem 1.1. Let X be a smooth vector field on a manifold M . Assume that the following 
three conditions hold: 

(i) There exists a vector fieldY G X°°(M) and two functions A,/i G C°°(M) such that [X,Y] = 
XX + fiY. 

[a) X admits an integrating factor f G (M) . 
[iii) The function f^ is an integrating factor ofY. 

Then the function 

H:=^ + divn{Y)-X (1.1) 

is a smooth first integral of the vector field X in the open set U := M\Zf, where Zj := {p G 
M:f{p) = 0}. 

Remark 1.1. The theorem holds true if X, Y and / are of class C^; in this case the first 
integral H is of class C^. On the other hand, if the zero set Zj has empty interior the first 
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integral H is defined in an open and dense subset of M. This is the case, e.g. if the integrating 
factor / is analytic (C^). 

Let us briefly explain the geometric meaning of the first integral derived in Theorem 1.1, cf. 
equation (1.1). We are very grateful to one of the referees for pointing out to us the argument 
that we sketch below. It is obvious that the function g := Hf is a second integrating factor of 
the vector field X; conversely, if X admits two integrating factors / and g it is ready to check 
that the function I := g / f is a first integral of X, but it does not need to be the same as the first 
integral H computed in Theorem 1.1 (see e.g. Example 2.1 in Section 2). The point is that under 
the assumptions of Theorem 1.1 it is possible to compute a second integrating factor g using 
the vector field Y and the integrating factor /, such that the function g/f is precisely the first 
integral H given by equation (1.1). For example, in the particular case that y is a symmetry 
of X (i.e. [X, y] = 0), then g := Lyf is a second integrating factor of X, thus implying that 
H := ^^y^ is a first integral of the vector field X. In the general case, one can look for a second 
integrating factor of the form g := Lyf+g, and Theorem 1.1 provides some sufficient conditions 
to compute the unknown function g. This is the main idea behind the proof of Theorem 1.1, 
which is given in Section 2. 

Let us also observe that the first integral H defined in (1.1) could be trivial (a constant), in 
which case our theorem does not provide any new information about the vector field X . The 
fact that this case can indeed happen (even when X has non-trivial first integrals) is illustrated 
in Section 2, thus showing that our way of producing first integrals is by no means the most 
general one. Finally, in Section 3 we give some applications and examples of our main result 
and show how it generalizes some previous theorems in the literature. 

2 Proof of Theorem 1.1 and an example 

First we shall prove Theorem 1.1. 

Proof of Theorem 1.1. Our method of proof, which is coordinate- free, makes use of standard 
tools from differential geometry, the reader who is not familiar with this differential geometric 
setting can consult e.g. [1, Chapter 2]. We first recall the following identity: 

^[X,Y] ~ ^xLy - LyLx- (2.1) 

Using the properties of the Lie derivative and assumptions (ii) and (Hi) in Theorem 1.1, it 
is easy to check that 



On the other hand, the identity (2.1), assumption (ii) and the fact that the Lie derivative 




of the Lie derivative implies that 



L[x,Y]if^) = Lxxim + L^Yim = Lx{X)fn. 



(2.4) 
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which imphes that H is a first integral of the vector field X. Observe that H is well defined 
(and C°°) in the set U := M\Zf, where Zf := {p £ M : f{p) = 0}. This completes the proof of 



Remark 2.1. Arguing exactly as in the proof of Theorem 1.1 it is obtained that the function 



is a second integrating factor of the vector field X. This implies that H := g/f is a first integral 
of X, thus providing an alternative approach to the proof of Theorem 1.1 in the spirit of the 
explanation given at the end of Section 1: the first integral H follows from the existence of an 
appropriately chosen second integrating factor. 

We finish this section by observing that, in the case that the function H is identically constant, 
it provides a trivial first integral of X. The following example shows that this situation can 
indeed happen, although it is not the general case as we shall see in Section 3. 

Example 2.1. Let X be a volume-preserving vector field having a volume-preserving centrali- 
zer Y which is (almost everywhere) transverse to X, that is 



It is straightforward to check that assumptions {i)-{ui) in Theorem 1.1 are satisfied by the 
vector field X. On the other hand, the function H is identically zero, thus providing a trivial 
first integral of X. When the manifold M is 3-dimensional, it can be shown [9, 15] that 

d{ixiY^) = 0, 

thus concluding, if the first De Rham cohomology group of M is trivial, that is H^{M;M.) = 0, 
that there exists a function / G C°°(M) such that dl = ixiy^- This function / is non-trivial 
because the vector field Y is transverse to X, and since Lx{I) = ixdl = ixixiy^ = 0, it is 
a first integral of X. Accordingly, the vector field X has a global first integral which is not 
produced by the mechanism involved in Theorem 1.1. 

Remark 2.2. Note that we have not assumed that the vector field Y be (almost everywhere) 
transverse to X in the statement of Theorem 1.1. If y is proportional to X, i.e. Y = hX for some 
function h, then we have that A = Lx{h), and being divQ(X) = — easily follows that 
H = 0. We conclude that, in order to obtain a non-trivial first integral H with Theorem 1.1, 
the vector field Y must be transverse to X. 

3 Applications and examples 

A particularly useful application of Theorem 1 . 1 is when the vector field y is a normalizer of X 
(a particular case of assumption (i) in the theorem) and X admits an integrating factor. In 
this case, a well known property of the normalizers of a vector field allows to obtain many first 
integrals (not necessarily independent). The following corollary is a generalization of a theo- 
rem in [6, 7] which states that divn(y) is a first integral of X provided that [-'^, y] = and 



the theorem. 



g:=LYf + fdwn{Y)-Xf 



[X,Y] = 0, 



divn(X) = dwn{Y) = 



rank(X, Y) = 2 a.e. in M. 



divniX) = 0. 
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Corollary 3.1. Assume that a vector field X G X°°(M) has a normalizer Y G X°°{M), 
i.e. [X,Y] = XX for some function A G C°°(M). Then, if X admits an integrating fac- 
tor f, the function H defined in (1.1) is a first integral of X . Moreover, the function '■= 
Ly{Ly{- ■ ■ {Ly{H)))) is also a first integral of X , where the Lie derivative Ly is applied k > 1 
times on H . 

Proof. It is obvious that the hypotheses of Theorem 1.1 are fulfilled with fi = 0, and therefore 
the function H is a first integral of X. If we prove that the function Hi := Ly{H) is a first 
integral of X it easily follows by induction that any function is also a first integral of X, so 
let us prove that this is the case for Hi. Indeed, 

Lx{Hi) = Lx{Ly{H)) = L^x,y]{H) + Ly{Lx{H)) = \Lx{H) = 0, 

where we have used the identity (2.1) and that Lx{H) = 0, thus proving the claim. ■ 

Notice that at most n — 1 independent first integrals can be obtained with this procedure, 
so most of the functions Hk will not add any new information. On the other hand, it is relevant 
to remark that the existence of an integrating factor for a vector field is not enough in general 
to guarantee the existence of a first integral, so it is necessary to assume further hypotheses, as 
in the statement of Corollary 3.1. 

In the following corollary we give a sufficient condition in order that a vector field be com- 
pletely integrable, i.e. to admit n — 1 independent first integrals (we say that two functions are 
independent if they are functionally independent in an open and dense subset of M). Com- 
pletely integrable vector fields were studied in [17] concerning their connection with the period 
function and the existence of normalizers. Let us observe that this result is quite different from 
a theorem recently proved by Kozlov [12]: on the one hand, we do not need to assume that 
the set of vector fields {X,Yi, . . . ,Yn-2} generates a (solvable) Lie algebra and on the other 
hand the first integrals that we obtain are explicit (i.e. divQ{Yi)) and there is no need of using 
quadratures to compute them. 

Corollary 3.2. Let X £ X°°(M) be a vector field whose flow preserves the volume form 0, and 
admits n — 2 normalizers Yi, . . . , Yn-2 so that [X, Yi] = CiX with CjGM, i=l,...,n — 2. Then, 
if the functions div^lYi) are independent, the vector field X is completely integrable. 

Proof. By assumption, divn(-'^) = 0, and hence Corollary 3.1 can be applied with / = 1 and 
A = Ci to imply the existence of n — 2 first integrals Hi = divQ(li). Under the assumption that 
the functions Hi are independent, it is standard that the pull-back of the vector field X, which 
is volume-preserving, to the regular level sets (of dimension 2) of the map (Hi, . . . , Hn-2) '■ 
M W--^ admits an integrating factor with respect to the induced volume form, and hence it 
is integrable [8, Chapter 2.11]. This defines a first integral Hn-i in any contractible domain of 
the regular level sets, thus proving the claim. ■ 

Remark 3.1. According to Corollary 3.2, the vector field X has n — 2 first integrals given 
by Hi = divfi{Y) (i = ^, ■ ■ ■ ,n — 2), and hence we can apply [2, Chapter 4.4, Theorem 13] 
to conclude that on each regular level set of the first integrals the field can be integrated by 
quadratures. Moreover, assuming that the regular level set is connected and compact and X 
is non-vanishing, it follows that it is diffeomorphic to a torus and X is orbitally conjugated to 
a linear field of (rational or irrational) frequency on the level set. Observe that, in order to 
use [2, Chapter 4.4, Theorem 13], it is necessary to have n — 2 first integrals, and Corollary 3.2 
precisely gives a sufficient condition for this. 

Theorem 1.1 can be applied to compute first integrals of Hamiltonian vector fields with 
symmetries (note that a Hamiltonian vector field is divergence-free with respect to some volume 
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form fi). In the following example, which is inspired by a paper of Hojman [10], we apply it to 
the motion of a particle in a spherically symmetric homogeneous potential to obtain the well 
known conservation laws of the energy and the angular momentum. We find this result rather 
surprising, as the method of proof is completely different from the classical one, thus showing 
the existence of a mechanism different from Noether's theorem to produce first integrals in 
Hamiltonian mechanics. 

Example 3.1. Let V{u) be a real valued function which is homogeneous of degree a G M\{0}, 
that is V{u) := u". Endowing M^" with Cartesian coordinates {xi, . . . ,Xn,Pi, ■ ■ ■ ,Pn), the 
Hamiltonian vector field which describes the dynamics of a particle under the action of the 
potential V{xf + • • • + x^) is given by 

X = pidx^ H ^Pndx„ - 2xiV'dp^ 2xnV'dp^, 

where V'{u) = avP^~^ . It is obvious that X is divergence- free with respect to the standard 
volume form = dx\ A • ■ ■ A dxn A dpi A ■ ■ ■ A dpn, that is div^X = 0. A straightforward 
computation shows that the vector field 

Yij := {xiPj - XjPi)[a'^{xidxi H h XnO^J + Pidp^ H V Pndp,,] 

verifies 

[X,Yij\ = (a^^ - l){xiPj - XjPi)X 

for all 1 < i < j < n. Applying Theorem 1.1 with / = 1 and Xij 
conclude that X has the following set of first integrals: 

Hij = divn(lij) - Xij = (n(l + a'^) + 2){xiPj - xjpi), 

for all 1 < i < J < n, which just expresses the well known conservation law of the angular 
momentum. Observe that the first integral Hij is trivial only when a = ^^r;^- Analogously, if 
we consider the vector field 

Y:= I^I{pI + ---+pI)+V{xI + --.+xI)'^ 

X [a''^{xidxi H \-Xnd:,„) +pidp^ H \- Pndp„] 

we get that 

[X, Y] = (a-i - 1) Q (p2 + . . . + p2^) + y (^2 ^ . . . ^ ^2)^ 

and hence Theorem 1.1 implies that X has the following first integral as well: 

H = (n(l + a-') + 3 - a'') Q (p? + • • • + pI) + V{xl + ■■■ + xl)^ . 

This first integral reflects the conservation of the energy, and is non-trivial unless a = 
Note that this example includes the n-dimensional harmonic oscillator (with a = 1) and the 
n-dimensional Kepler problem (with a = (2 — n)/2) for n > 2. 

We finish this letter by elaborating on the 2-dimensional case. This is a rather exceptional 
situation because the existence of a normalizer implies the existence of an integrating factor, 
and vice versa, and therefore Theorem 1.1 can be applied. In Example 3.2 below we show that 
the first integral H in (1.1) turns out to be trivial for the usual choice of integrating factor 
(normalizer) constructed from a given normalizer (an integrating factor). Nevertheless, the first 
integral H provided by Theorem 1.1 is not trivial in general for 2-dimensional vector fields, as 
shown in Example 3.3. 



= (a ^ - l){xiPj - XjPi) we 
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Example 3.2. Let X be a vector field on a 2- dimensional manifold M which admits a transverse 
normalizer Y , i.e. [X, Y] = \X and rank(X, Y) = 2 a.e. in M. It is immediate to check [5] that 
the function / := {ixiY^)~^ is an integrating factor of X. We now compute the function A in 
terms of Y and /. Defining the vector field W := fX and the 1-form oj := iw^, it follows that 
iiyiyO = 1 and using Cartan's formula for the Lie derivative we get 

L]y{iY^) = iwidiy^) + d{iwiY^) = divQ{Y)uj. (3-1) 
On the other hand, since the identity L^iy = i[w,Y] + iyLw holds, we derive that 

Lw{iYn)=(^X-^^'^u;, (3.2) 

where we have used that [IV, = — ^^^y^^ W and that Lw{^) = 0. Identifying equa- 
tions (3.1) and (3.2) and noticing that the 1-form uj is not trivial, we conclude that 

A = div.(y) + MZ). 

Therefore, the first integral H defined in (1.1) is identically zero for this choice of integrating 
factor / given a normalizer Y. 

Conversely, if X admits an integrating factor /, it is easy to check [5] that the vector field 
Y := [f'^ divQ{X))~^ Z is a normalizer of X, with Z the unique vector field which solves the 
equation iz^ = —df. Arguing exactly as before it is straightforward to check the normalizer 
condition [X, y] = XX, where 

A = divo(y) + ^, 

and therefore the first integral if of (1.1) is also trivial for this choice of normalizer Y given an 
integrating factor /. 

We observe that a non-trivial first integral of X can be obtained through the classical Lie's 
integration algorithm [16] for two dimensional vector fields with a transverse normalizer Y. 
Indeed, the 1-form 

UJ 



ix^Y^ 



is closed, and hence there exists a smooth function / defined in any contractible domain of the 
manifold M such that dl = uj. It is clear that this function / is a first integral of the vector 
field X. 

Example 3.3. Consider the following polynomial vector field in (we endow the plane with 
Cartesian coordinates (x, y) and the standard volume form dx A dy): 

^ = rr^ [-^y'^- + {i + x'+y'+ 4x^)9,] . 

It is easy to check that 

[X,Y] = XX, 
where the vector field Y is given by 

Y = i^ + ^' + y'?f [(1 + + / + 4x3)5. + Ay^dy] , 

(1 + x4 + y4 + 4^3)2 + igy6 ! 
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and the expression of the function A is omitted for the sake of simphcity. Moreover, the function 



is a non-trivial global first integral of X. 

For the sake of completeness, we finally compute the first integral I{x, y) of X obtained 
through the classical Lie's integration algorithm (see the last paragraph of Example 3.2). After 
a few computations we get that 



thus showing that, in general, our approach gives a first integral different from the one obtained 
through classical methods, although of course in the two-dimensional case both first integrals 
are functionally dependent. 
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